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比值 校正 法 在 射电 频谱 中 的 应 用 
张 燕 坤 ”"， 董 亮 ， 李 明 悦 ， 何 乐 生 ， 汪 敏 ” 
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摘要 : er co id ei 
差 ， 难 以 获得 分 子 谱 线 准确 的 频谱 信息 。 针 对 射电 频谱 上 离散 数字 功率 谱 分 析 中 带 来 的 频率 偏 
移 问 题 进行 比值 校正 法 研究 和 处 理 ， 该 方法 是 频谱 矫正 方法 的 一 种 ， 特 点 在 于 算法 简单 ， 方 便 
快捷 。 利 用 射频 信号 相同 、 分 辨 率 不 同和 不 同 射频 信号 、 不 同 分 辨 率 情 况 下 ， 归 一 化 频率 偏 移 
量 的 变化 确定 该 方法 的 校正 效果 , 最终 实现 了 ” 0MHz, 141. 8MHz 和 270. 8MHz 的 频谱 校正 ， 
提高 了 相应 的 频谱 分 辨 率 可 达 7. $SHz/channel。 通 过 实验 结果 对 比 ， 确 定 了 测试 系统 的 最 佳 频 率 
分 辨认 为 50Hz/Channel 或 200Hz/channel。 
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谱 线 在 研究 天 体 的 物理 或 化 学 环境 、 演 化 历程 以 及 天 体 分 类 等 方面 具有 重大 作用 ， 是 天 文 
学 研究 的 一 个 工具 。 与 各 种 辐射 机 制 预言 的 功率 谱 相 似 ， 谱 线 研 究 涉及 一 个 频率 、 多 个 频率 或 
者 连续 谱 ， 可 以 是 吸收 线 ， 也 可 以 是 发 射线 。 它 们 来 自 于 不 同 层级 的 天 体 ， 比 如 恒星 、 星 云 、 
星系 、 本 星系 群 或 者 整个 宇宙 。 通 过 多 波段 观测 可 以 对 这 些 天 体形 成 一 个 系统 的 认识 。 在 众多 
的 谱 线 中 , Hi21 厘 米 谱 线 是 验证 大 爆炸 理论 ， 揭 示 宇宙 早期 再 电离 阶段 的 重要 探 针 ”; 光谱 学 中 
的 各 类 线 系 、 复 合 线 可 以 用 来 研究 天 体 的 温度 、 速 度 等 物理 量 ， 人 尘埃 辐射 和 分 } 子 谱 线 是 研究 恒 
星 形成 、 星 际 介质 的 重要 工具 或 手段 。 

各 类 谱 线 具 有 不 同 的 波长 ， 其 中 很 多 在 射电 波段 内 〈0. 3mm 一 20n) ， 可 以 通过 射电 望远镜 
或 天 线 阵 列 对 它们 进行 探测 。 为 了 尽 可 能 探测 到 较 多 的 谱 线 ， 除 了 需要 好 的 射电 天 文 观 测 环境 
外 ， 还 要 望远镜 的 灵敏 度 和 分 辨 率 足 够 好 。 因 此 ， 低 噪声 放大 器 、 宽 于 接 收 机 、 大 口径 天 线 、 
多 波束 或 焦 面 阵 《〈 大 视 场 ) 等 关键 设备 的 研发 成 为 国际 和 国内 的 重大 课题 。 同 样 ， 多 通道 数字 
终端 、 高 速 数据 处 理 和 存储 设备 研发 ， 甚 至 观测 技术 或 方案 等 也 很 重要 。 

除了 空间 分 辩 率 以 及 时 间 分 辩 率 的 要 求 ， 分 子 天 体 物理 和 谱 线 观测 对 终端 设备 的 频率 分 辨 
率 也 提出 了 很 高 的 要 求 。 在 星际 介质 和 恒星 形成 研究 工作 中 ， 天 体 的 运动 速度 或 速度 梯度 场 

是 很 重要 的 物理 量 。 它 在 天 体 测 量 学 、 星 系 动力 学 中 定义 的 本 地 静止 标准 中 表现 为 红 移 或 蓝 移 ， 
即 多 普 勒 运动 这 个 速度 会 使 探测 到 谱 线 频率 相对 于 地 球 上 实验 室 测 得 的 频率 有 一 定 漂移 ， 此 
时 终端 较 高 的 频率 分 辨 率 会 使 这 样 的 漂移 被 精确 测定 ， 从 而 提高 速度 测量 的 可 信和 度 。 另 外 ， 展 
宽 也 就 可 以 定 得 比较 准确 。 如 图 1， 从 H(96) B 的 离散 谱 形 能 够 得 到 天 体 的 速度 、 温 度 等 物理 信 
息 ， 表 明 频 谱 分 辩 率 越 高 ， 相 关 研 究 越 能 反映 真实 的 天 体 物理 信息 。 
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Velocity (km/s | 
图 1 天 马 65 米 望 远 镜 观 测 到 的 H(96) B 频谱 。 图 像 上 沿 坐 标 轴 显 示 的 是 经 过 多 普 勒 修正 后 的 频率 ， 以 MHz 为 单位 ; 纵 轴 表 示 天 线 温 


度 ， 以 K 为 单位 。 


Fig. 1 The H(96) B spectrum observed by Tianma 65 meters telescope. The upper horizontal axis shows the frequency in unit of 
MHz after Doppler correction. The vertical axis represents the antenna temperature in unit of K. 


实际 上 ， 观 测 设备 和 观测 方法 会 在 频谱 上 造成 一 定 的 频率 漂移 和 速度 展 宽 ， 从 而 引入 系统 
误差 ， 因 此 还 需要 频谱 校正 。 尽 管 该 技术 在 工程 领域 得 到 广泛 应 用 ， 但 是 天 文 领域 的 应 用 仍然 


欠缺 。 在 频谱 分 析 时 ， 一 般 采 用 拟 合 方法 检验 通道 稳定 


性 和 通道 间隔 


致 性 ， 很 少 用 到 频谱 校 


正 的 方法 ,也 没有 校正 偏 移 量 。 这 对 于 频率 分 辨 率 要 求 不 高 (6~10 一 10 Hz ) 的 谱 线 观测 是 合适 
性 也 相应 提高 ， 此 时 简单 的 拟 


的 ， 但 是 ， 随 着 天 文 观 测 要 求 的 提高 ， 接 收 机 、 终 端 精 
合 不 再 适应 ， 频 谱 校正 就 显得 尤为 重要 。 


度 和 稳定 4 


频谱 校正 是 信号 处 理 领 域 重 要 的 研究 课题 。 
行 计算 ， 因 此 ， 需 要 对 时 域 信号 进行 


由 于 计 


算 机 只 能 对 有 限 多 个 样本 和 采样 点 数 进 


断 和 数字 离散 化 ， 这 不 可 避免 地 造成 了 能 量 泄漏 和 栅栏 


效应 , 使 得 离散 频谱 的 幅 值 、 相 位 和 频率 产生 较 大 的 误差 *。 在 对 单 频率 谐 波 信号 加 窗 处 理 后 ， 


幅 值 误 差 最 高 可 达 36. 4%， 相 位 误差 可 达 土 90”， 频 率 误 差 可 达 圭 0.5 个 频率 分 辨 率 门 。 


国内 外 对 离散 频谱 校正 的 方法 有 很 多 ， 其 中 


目前 ， 


比值 校正 法 较为 常用 。 根 据 原 时 域 信号 


组 成 结构 的 不 同 ,主要 分 为 两 类 :一 类 频谱 校正 方法 主要 针对 


于 射电 天 文 的 谱 线 信号 主要 是 单 频 信 号 ， 即 使 是 能 级 分 


单 频 信号 , 另 一 类 针对 密集 频谱 。 由 


裂 所 造成 的 各 点 频 信号 间隔 也 较 大 ， 因 
此 ， 本 文 只 讨论 单 频 信 和 号 的 频谱 校正 。 对 于 单 频 信 号 ， 国 内 外 有 5 种 对 功率 谱 校 正 的 方法 : 比 


值 校正 法 、 能 量 重 心 校正 法 、FFT+FT 连续 细 化 分 析 傅 上 


FT 法 中。 本 文 的 实验 模拟 了 射电 天 文 观测 软 硬 件 系统 ， 
高 的 比值 校正 法 中 ， 对 经 过 LabView 程序 、 捷 变 收发 器 


1 比值 校正 法 


在 谱 线 研究 中 ， 频 率 、 流 量 和 偏振 是 比较 重要 的 观 


E 叶 变换 法 、 相 位 差 法 、 相 位 差 法 + 单 点 


并 采用 算法 简单 、 运 算 速 度 快 、 精 度 较 


的 噪声 源 信 号 进行 频率 定 标 。 


测量 。E 


分 观测 中 会 间接 影响 流量 乃至 偏振 的 测量 与 定 标 ， 


用 。 


于 频率 稳定 性 和 测量 偏差 在 积 


因此 , 仅 介绍 比值 校正 法 在 频率 定 标 方面 的 应 


利用 频率 归 一 化 后 差 值 为 1 的 主办 峰 项 附近 两 条 谱 线 的 窗 谱 函 数 比 值 , 建立 一 个 以 归 一 化 校 


正 频率 为 变量 的 方程 ， 
校正 法 外 。 解 出 归 一 化 校正 频率 的 方法 有 多 种 ， 


解 出 归 一 化 校正 频率 ， 进 而 进行 


直接 时 


频率 、 幅 值 和 相位 校正 的 方法 称 为 比值 


公式 的 方法 称 为 比值 公式 法 ， 利 用 连 


代 求 解 的 方法 称 为 比值 迭代 公式 法 ， 用 搜索 求解 的 方法 称 为 比值 峰 搜 索 法 。 


假设 归 一 化 窗 函 数 的 频谱 模 函 数 表达 式 为 W10)， 


于 主办 内 


幅 值 为 了 


E， 此 时 频谱 模 函 数 与 


归 一 化 窗 谱 函 数 完全 相等 ，Wi( 几 对称 与 y 轴 ， 即 主 匆 中 心 是 坐标 原点 。 
如 果 周 期 信号 的 频率 恰好 对 正 某 一 发 射 ( 吸 收 》 谱 线 频率 ， 那 么 计算 得 到 的 频率 、 幅 值 和 


相位 都 是 准确 的 。 更 一 般 的 情况 是 ， 信 号 频率 没落 在 通 


道中 心 位 置 ， 因 


此 ， 峰 值 谱 反 映 的 频率 
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和 幅 值 都 不 准确 ， 相 位 误差 也 较 大 。 主 办 中 心 并 不 在 坐标 原点 ， 而 是 有 一 定 偏 移 W(-0.5<VA < 


0.5)。 


在 导电 天 文 的 高 精度 实测 中 ， 需 要 精确 测量 得 到 谱 线 的 频谱 信息 ， 因 而 频谱 的 偏 移 要 尽 可 


能 消除 。 


设 和 表示 谱 线 发 射频 率 ， 大 表示 频谱 仪 接收 的 频率 ， 则 在 设备 恨 好 的 情况 下 ， 地 面 测试 实验 
得 到 的 结果 应 是 =fr。 但 是 ， 随 着 模 数 转 换 技术 的 发 展 和 离散 数字 信号 的 应 用 ， 连 续 的 频谱 被 
离散 化 。 由 于 谱 线 发 射频 率 处 在 特定 通道 K( 幅 值 谱 谱 线 号 ) 内 ， 但 不 一 定 处 在 该 通道 的 中 心 
位 置 处 ， 因 此 ho。 加 上 时 域 信号 被 截断 的 影响 ， 与 1 之 间 的 偏差 较 大 。 经 过 频谱 校正 ， 可 


以 在 较 高 精度 (Hz 级 ) 下 ， 可 以 使 得 fh .==f。 

实际 频谱 中 ， 假 设 采样 点 数 为 2V， 采 样 频率 为 2 大 ， 那 么 通道 数 为 W， 采 检 
KE[1,VW2-1]5。 若 通道 K 的 幅 值 为 下 ， 其 两 侧 K 1、K+1 通 道 的 幅 值 分 别 为 双修 
所 对 应 的 频率 为 


[六 


f=K - (1) 
当 > 时 ， 如 图 2， 有 
dd 
C= (L E[-0.5,0]) (2) 
了 
当 Za<Za 时 ， 如 图 3， 有 
过 
2 = 一 (LE[0,0.5]) (3) 
Yes 
对 于 和 矩形 窗 而 言 ， 其 频率 校正 量 为 
1 
Yr > 六 时， We 
1 (4) 
Ye > Yr 时 ,Vf = 一 二 
1 十 一 
(0 
校正 得 到 的 校正 频率 5 为 
万 = (K+Vf') (5) 
N 


EF 带 宽 为 AR， 且 
有 1。 那 么 有 通道 


如 图 4， 以 上 校正 算法 适用 于 起 始 频率 为 0Hz 的 离散 频谱 ， 但 实际 应 用 中 ， 该 起 始 频率 一 般 不 
为 0Hz， 因 此 ， 还 需要 针对 特定 领域 进一步 计算 。 在 射电 频谱 中 ， 真 实 的 频率 为 


f=fiotfo ec (6) 


其 中 ，fio 表 示 本 地 振荡 占 的 频率 。 千 有 多 级 本 振 ， 还 应 把 各 级 本 振 频 率 求 和 ， 才 能 得 到 fio。 


202004.00029v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


Psis 
ideal 
S case 1 
富 EE 
om 
Pp 
5 1| 
目 
© EL 
0 
0 , | 
0 flueal F 
Frequency (Hz) 
图 2 VV 时 对 应 的 离散 信号 示意 图 
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图 4 考虑 本 地 振荡 频率 后 的 校正 示意 图 


Fig.4 _ Calibration schematic diagram after considering the local oscillation frequency 


2 实验 设置 


射电 天 文 观测 中 ， 天 文 谱 线 信号 在 L 波段 及 其 以 上 频率 (>1GHz)， 因 此 需要 降 频 ， 以 降低 
数字 后 端 采样 的 压力 。 观 测 源 的 信号 经 过 天 线 、 接 收 机 后 ， 其 频率 被 调制 到 几 十 到 几 王 兆赫 效 ， 
随后 将 这 样 的 中 频 信 号 送 至 终端 进行 采样 和 校准 。 

为 了 模拟 出 这 样 的 一 个 中 频 信号 ， 采 用 标准 信号 源 并 设置 输出 频率 为 


CNINa IV1 人 | F 其 


刻 =141,800,000.00Hz。 同 时 ， 使 用 捷 变 收发 器 NI2901 作为 模 数 转换 和 数据 采样 设备 ， 并 编写 
LabView 程序 设置 终端 参数 和 控制 运行 过 程 。 男 外 ， 通 过 程序 中 的 功率 谱 计 算 ， 在 输入 带宽 内 
划分 出 成 驻 个 数字 通道 从 而 将 一 定 的 信号 带宽 细 化 ， 这 样 就 能 提高 频率 分 辨 率 ， 满 足 对 频谱 
分 辩 率 的 不 同 需求 。 
实验 设置 了 不 同 的 频率 分 辨 率 ， 以 观察 比值 校正 法 在 不 同 天 文 需求 下 的 表现 。 虽 然 分 辩 率 
越 高 对 天 文 观 测 越 有 利 ， 但 是 需要 与 当前 技术 发 展 水 平 相 适 应 ， 因此 参数 配置 要 合理 » 具体 情 
况 见 表 1。 从 表 1 中 载波 输入 和 实际 载波 部 分 可 以 看 出 ， 设 定 的 参数 与 设备 响应 间 有 一 定 的 差 
值 ， 这 给 数据 分 析 带 来 不 便 。 同 时 ， 为 了 使 实验 具有 一 定 的 统计 意义 ， 对 每 组 设置 分 别 测试 5 
次 。 


表 1 不 同 频 率 分 辩 率 下 的 参数 设置 


Table 1 Parameter Settings at different pt resolutions 


Freq. Resol.” |0 yate Sig. BW Act. BW Finput” Factual" Sampling Cumulative 
(Hz/channel) Eo (Hz) (Hz) (Hz) Wie points frequency 
站 00. 0OK 22. 00K 22. 00 141769500 141. 769 16384 1000 
15 00. 00K 22. 00K 22. 00 41769500 41.769 8192 1000 
20 00. 00K 81. 92K 81. 92K 141769500 41.770 4096 1000 
30 00. 00OK 22. 00K 22. 00 141769500 41.770 4096 1000 
50 00. 00K 22. 00K 22. 00 41769500 41.770 2440 1000 
100 00. 00K 22. 00K 22. 00 41769500 41.770 1220 1000 
200 00. 00K 22. 00K 22. 00 141769500 41.770 610 20000 
500 00. 0OK 22. 00K 22. 00 141769500 141.770 244 20000 
1000 00. 00K 22. 00K 22. 00 41769500 41.770 122 20000 


Frequency resolution; Signal band-width; Actual band-width; Input carrier frequency; Actual 


carrier frequency. 


>q 


寸 于 特定 观测 波段 ， 由 于 观测 源 和 探测 的 分 子 不 同 ， 其 输出 的 中 频 信 和 号 也 不 同 。 考 虑 到 40 
米 射 电 望 远 镜 分 子 谱 线 终端 的 中 频 输 入 带宽 是 1GHz (升级 改造 后 为 2GHz )， 谱 线 信 号 则 可 能 
出 现在 带 通 的 几 十 或 几 百 MHz 的 位 置 ， 因 此 还 设置 了 80.0MHz 和 270.8MHz 的 实验 组 ， 以 丰 
究 比 值 校正 法 在 不 同 输出 信和 号、 不同 频率 分 辨 紊 下 的 表现 。 与 表 1 参数 不 同 的 是 ， 当 射频 信号 
为 80MHz 时 , 载波 输入 为 79.9795MHz, 实际 载波 为 79.9795MHz; 当 射 频 信 号 为 270.8MHz 时 ， 
载波 输入 为 270.7695MHz， 实 际 载 波 为 270.77MHz。 


3 结果 分 析 


过 多 次 测试 和 数据 分 析 ， 发 现 和 解决 了 理论 和 技术 上 的 缺陷 ， 并 有 了 一 些 结果 。 认 真 分 
析 这 些 结果 ， 有 助 于 进一步 研究 频谱 校正 理论 ,并 把 相关 的 方法 应 用 到 实际 观测 中 。 本 文 针对 不 
同 目的 设置 了 不 同 的 实验 组 ， 下 面 分 别 叙述 。 

首先 ， 归 一 化 频率 偏 移 量 VF 在 射频 信号 相同 ， 频 率 分 辨 率 不 同 的 情况 下 表现 有 一 定 规律 。 
在 表 2 和 图 5 中 ，0~4 分 别 表示 5 次 测试 的 序号 。 从 图 中 可 以 看 出 ， 当 频率 分 辩 率 太 高 时 ， 各 
次 测试 的 值 不 尽 相 同 ， 而 在 较 低 频率 分 辨 率 时 各 次 测试 一 致 性 较 好 ， 这 说 明 系 统 在 较 高 频谱 分 
状 率 要 求 下 不 够 稳定 ,并 且 WF 大 于 0 的 散 点 较 多 ， 即 多 数 情况 是 Yi.1<Yiw1， 这 与 表 3 和 图 6 相 吻 
合 。 另 外 ， 5 和 图 6 还 表明 : 不 论 分 辨 紊 的 高 低 ， 归 一 化 频率 偏 移 量 总 在 [-0.5,0.5] 之 间 。 因 
此 ， 比 值 校正 法 在 较 高 频率 分 辨 率 下 的 校正 效果 更 好 。 对 于 本 实验 所 用 的 设备 ， 实 际 观测 时 的 
频率 分 辩 率 可 在 50Hz/channel，100Hz/channel 或 200Hz/channel 间 选 择 。 
其 次 ， 频 率 分 辨 率 相同 ， 但 是 射频 信号 不 同 ， 则 归 一 化 频率 偏 移 量 VA 也 不 同 。 从 表 4 和 图 
7 可 以 看 出 ， 射 频 信号 分 别 为 80MHz，141.8MHz 和 270MHz 时 ， 归 一 化 频率 偏 移 量 有 较 大 差 
别 , 并 且 频 率 分 辨 率 较 高 和 较 低 时 有 很 大 不 同 。 从 图 中 可 以 看 出 : 当 频 率 分 辨 紊 为 50Hz/channel 


时 
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或 200Hz/channel 时 ，3 个 信号 的 归 一 化 频率 偏 移 量 YP es [-0.2，0.2]。 结 合 图 5、 图 6， 可 知 此 时 


系统 最 稳定 。 


以 上 结果 是 可 以 理解 的 : 由 于 比值 校正 法 只 使 月 
很 高 的 依赖 性 。 如 果 在 高 频率 分 辩 率 情况 下 观测 ， 


法 是 最 好 、 最 快 的 校正 方法 。 


了 3 个 临近 点 ， 


有 分析 


归 一 化 频率 偏 移 


表 2 射频 为 141.8 MHz， 频 率 分 辩 率 不 同时 的 多 次 测试 数据 


内 此 ,对 频率 分 辩 率 显示 出 
旦 ， 那 么 比值 校正 


T 


Table 2 Results with radio frequency of 141. 8MHz and different frequency resolution 


Index of 


Fred. 


resol. 


Index of 


Fred. 


resol. 


Vr VA 
exper iments (Hz/channel) experiments (Hz/channel) 
1000. 0000000000 0. 3423212243 1000. 0000000000 0. 0651720930 
100. 0000000000 0. 0142455346 100. 0000000000 0. 0144940334 
7. 4462890625 -0. 0416315433 7. 4462890625 0. 0244703661 
500. 0000000000 0. 3508046810 500. 0000000000 0. 3487047748 
0 50. 0000000000 0. 0387933055 3 50. 0000000000 0. 0570070349 
29. 7851562500 0. 2493147110 29. 7851562500 0. 3044267795 
200. 0000000000 0. 0047153864 200. 0000000000 0. 0046528023 
14. 8925781250 0. 3275106419 14. 8925781250 0. 4384830843 
20. 0000000000 -0. 0530305312 20. 0000000000 -0. 0706707358 
1000. 0000000000 0. 0655970348 1000. 0000000000 0. 0651342265 
100. 0000000000 0. 0155808114 100. 0000000000 0. 0185485020 
7. 4462890625 0. 0041361850 7. 4462890625 0. 3064197324 
500. 0000000000 0. 3494843586 500. 0000000000 0. 3423212243 
1 50. 0000000000 0. 0533924633 4 50. 0000000000 0. 0598957511 
29. 7851562500 0. 2645479331 29. 7851562500 0. 4164752968 
200. 0000000000 0. 0047438926 200. 0000000000 0. 0045723524 
14. 8925781250 0. 2875929373 14. 8925781250 -0. 4572477817 
20. 0000000000 -0. 0492047472 20. 0000000000 -0. 0470696696 
1000. 0000000000 0. 0653673598 
100. 0000000000 0. 0149879265 
7. 4462890625 0. 1447826447 
500. 0000000000 0. 3478842095 
2 50. 0000000000 0. 0541893805 


29. 7851562500 
200. 0000000000 
14. 8925781250 
20. 0000000000 


-0. 4956424680 


0. 0047095546 
0. 3512400507 
-0. 0473347745 
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表 3 射频 为 141. 8 MHz， 频 率 分 辨 率 不 同时 VA 的 均值 与 标准 差 
Table 3 The mean and standard deviation data withradio frequency of 141.8MHz and 
different frequency resolution 


Freq. resol. 


(Hz/channel) Mean of Wf Std. of WF 
100. 0000000000 0. 0155713616 0. 0015570680 
7. 4462890625 0. 0876354770 0. 1255422437 
1000. 0000000000 0. 1207183877 0. 1108015405 
200. 0000000000 0. 0046787977 0. 0000608743 
14. 8925781250 0. 1895157865 0. 3271404205 
50. 0000000000 0. 0526555871 0. 0072994991 
500. 0000000000 0. 3478398496 0. 0029223252 
20. 0000000000 -0. 0534620917 0. 0088640367 
29. 7851562500 0. 1478244505 0. 3270082220 


fres VS Vr 


0.41 


0.21 


0.0 ] 


二 人 2 


—0.4 4 


Te 


0 @r f=141.800000 MHz 
1@rf=141.800000 MHz 
2 @r f=141.800000 MHz 
3 @rf=141.800000 MHz 
4 @rf=141.800000 MHz 


T 
0 50 


T 
100 


T T T T T T 
150 200 250 300 350 400 450 


fres (Hz) 


T T T T T T T 
500 550 600 650 700 750 800 


图 5 射频 为 141. 8MHz， 频 率 分 辩 率 不 同时 的 多 次 测试 结果 


Fig.5 Multiple test results with radiofrequency of 141.8MHz and different frequency resolution 


T T T T 
850 900 950 1000 
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fies VS vf} 


一 100.000000 Hz/Channel @ r f=141.800000 MHz 
1.00 | 一 7.446289 Hz/Channel @ r f=141.800000 MHz 
一 上- 1000.000000 Hz/Channel @ r f=141.800000 MHz 
一 上- 200.000000 Hz/Channel @ r f=141.800000 MHz 
一 14.892578 Hz/Channel @ r f=141.800000 MHz 
一 50.000000 Hz/Channel @r f=141.800000 MHz 
一 全 500.000000 Hz/Channel @ r f=141.800000 MHz 
一 上 20.000000 Hz/Channel @ rf=141.800000 MHz 
0.50] 一 29.785156 Hz/Channel @ r f=141.800000 MHz 


0.75 1 


202004.00029v1 


chinaXiv 


0.25 1 


Vfl (Hz) 


0.00 1 说 * 


-0.25 ] 


-0.50 ] 
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T 
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图 6 射频 为 141. 8MHz， 频 率 分 辨 率 不 同时 的 分 析 结 果 


Fig.6 Analysis results of radio frequency about 141.8MHz with different frequency resolutions 


表 4 V1P 在 不 同 射频 ， 不 同 频率 分 辨 率 情况 下 的 值 
Table 4 The value ofVPwith different radio frequency and frequency resolution 


T T T T T T T T T T T T T T T T T T T 
0 50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 


RF Freq. Resol. vA RF Freq. Resol. ve 
(MHz) (Hz/channel) (MHz) (Hz/channel) 
1000. 0000000000 0. 2054970208 1000. 0000000000 0. 3423212243 
100. 0000000000 0. 2095411283 100. 0000000000 0. 0142455346 
7. 4462890625 0. 3715679928 7. 4462890625 -0. 0416315433 
500. 0000000000 -0. 1878904187 500. 0000000000 0. 3508046810 
80.0 50. 0000000000 -0.1722237239 141.8 50. 0000000000 0. 0387933055 


29. 7851562500 
200. 0000000000 
14. 8925781250 
20. 0000000000 


-0. 1263683613 
-0. 1672976339 
-0. 1497675356 
-0. 0128196597 


29. 7851562500 
200. 0000000000 
14. 8925781250 
20. 0000000000 


0. 2493147110 
0. 0047153864 
0. 3275106419 
-0. 0530305312 


1000. 0000000000 
100. 0000000000 
7. 4462890625 
500. 0000000000 
270.8 50. 0000000000 
29. 7851562500 
200. 0000000000 
14. 8925781250 
20. 0000000000 


-0. 0020558727 
0. 4999769741 
0. 4074170232 
-0. 0098171494 
0.0112350002 
-0. 3560928468 
-0. 0933705315 
-0. 0408522077 
0. 0727707923 
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4 总 结 与 展望 


比值 校正 法 简单 、 实 用 ， 这 给 财 电 天 文 


T T T T T T T T T T T T T T T T T T 
150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 


fres (Hz) 


图 7 不 同 射频 信号 ， 不 同 频率 分 辨 率 下 的 归 一 化 频率 偏 移 量 


Fig.7 Normalized frequency offsets at different radio signals and different frequency resolutions 


度 等 研究 工作 才 有 可 信 度 。 具 体 的 意义 有 3 点 : 


监测 可 以 得 出 它 的 变化 规律 ; 
本 振 的 频率 稳定 性 和 准确 性 ; 


该 工作 对 确定 仪器 误差 、 


频谱 分 析 带 来 了 很 大 的 便利 。 本 文 基于 测试 频 点 已 
知 的 情况 ， 虽 然 与 实际 观测 有 一 定 差距 ， 但 是 对 于 未 来 观测 数据 质量 的 提供 是 一 件 具 有 开拓 性 
的 工作 。 当 系统 误差 和 外 界 干扰 都 能 够 很 好 地 避免 时 ， 天 体 的 动力 学 或 运动 学 条 件 、 温 度 和 密 
(1) 检验 信号 源 或 者 本 振 的 频率 稳定 性 ， 长 期 


(2) 频率 切换 观测 中 ， 需 要 不 断 改 变 本 振 频率 ， 因 此 , 需要 监测 


(3) 实际 天 文 观测 中 ， 一 般 不 需要 对 频率 稳定 性 定期 检测 ， 但 是 


检测 仪器 性 能 以 及 对 更 深入 的 频率 定 标 和 监测 都 有 参考 价值 。 


在 未 来 ， 天 线 观测 谱 线 信号 的 频谱 前 ， 本 文 解决 了 频谱 校正 的 问题 ， 
的 分 辩 率 设置 (50Hz/channel 或 200Hz/channel) 给 动态 校 1 


同时 也 指明 了 实测 时 


E 频 谱 的 工作 提供 了 参考 。 另 外 ， 天 


体 的 多 普 勒 也 给 频谱 带 来 了 频率 漂移 ， 这 将 是 下 一 步 研究 的 内 容 ， 可 以 和 选择 观测 源 和 谱 线 等 


工作 一 起 进行 。 
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The Application of Ratio Correction in Radio Spectrum 


Yankun Zhang”, Liang Dong’, Mingyue Li’, Lesheng He’, Min Wang” 
(CUniversity of Chinese Academic of Sciences, School of astronomy and space science, Beijing100049, China) 
CChinese Academy of Sciences, Yunnan Astronomical Observatory, Kunming 650216, China) 
GYunnan University, School of Information, Kunming 650504, China) 

Abstract: Molecular Spectral lines can be viewed on the radio Spectrum, but discrete digital 
power spectrum analysis algorithmswould introduce errors that make it difficult to obtain accurate 
information about spectral lines. In this paper, we studied the ratio correction method and carried 
it out to solve the problem of frequency offset in discrete digital power spectrum analysis. This 
method is a kind of spectrum correction method, which is characterized by simple, convenient and 
fast algorithm.The correction effect of this method is determined by the change of the normalized 
frequency offset under the condition of the same radio signal, different spectral resolution and 
different radio signal, different spectral resolution. We finally complete the correction in tests of 
setting the radio frequency to be 80.0 MHz, 141.8MHz and 270.8MHz, and improve the frequency 
resolution to 7.5 Hz/channel. By comparing the experimental results, the optimal frequency 
resolution of the sampling system is determined to be 50.0 Hz/channel and 200.0 Hz/channel. 

Key words:Discrete digital power spectrum analysis;Ratio correction method;Spectral 
lines;spectral resolution 
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